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An investigation of dynamical chiral symmetry breaking on the light front is made in the Nambu- 
Jona-Lasinio model with one flavor and A'^ colors. Analysis of the model suffers from extraordinary 
complexity due to the existence of a "fermionic constraint," i.e., a constraint equation for the bad 
spinor component. However, to solve this constraint is of special importance. In classical theory, 
we can exactly solve it and then explicitly check the property of "light-front chiral transformation." 
In quantum theory, we introduce a bilocal formulation to solve the fermionic constraint by the 
1/N expansion. Systematic 1/N expansion of the fermion bilocal operator is realized by the boson 
expansion method. The leading (bilocal) fermionic constraint becomes a gap equation for a chiral 
condensate and thus if we choose a nontrivial solution of the gap equation, we are in the broken phase. 
As a result of the nonzero chiral condensate, we find unusual chiral transformation of fields and 
nonvanishing of the light-front chiral charge. A leading order eigenvalue equation for a single bosonic 
state is equivalent to a leading order fermion-antifermion bound-state equation. We analytically 
solve it for scalar and pseudoscalar mesons and obtain their light-cone wavefunctions and masses. 
AH of the results are entirely consistent with those of our previous analysis on the chiral Yukawa 
model. 

PACS number(s): ll.30.Rd, ll.30.Qc, 11.15.Pg, 12.40.-y 



I. INTRODUCTION 



Our expectation that light-front (LF) formahsm enables us to relate QCD directly to the constituent quark model 
at field-theoretic level Q, seriously requires a full understanding of dynamical chiral symmetry breaking (D^SB) on 
the LF. Central to this issue is, first of all, a well-known problem of how to reconcile a LF "trivial" vacuum with a 
chirally broken vacuum having a nonzero fermion condensate. The secondary problem is to determine the property of 
"LF chiral transformation" which is defined differently from the usual one. The most surprising fact of the LF chiral 
transformation is that it is an exact symmetry even for a massive free fermion ||^. 

In the present paper, we discuss this issue within the Nambu- Jona-Lasinio (NJL) model |^] which is a typical 
example of D^SB. Previously we consider the same problem from different point of view Q|. Our interest was in 
describing D^SB of the NJL model, but we actually took a roundabout way in order to apply an idea which works 
well for spontaneous symmetry breaking of a scalar model, to a fermionic theory. We know that the longitudinal zero 
modes of scalar fields are responsible for describing spontaneous symmetry breaking on the LF. Indeed, it is achieved by 
solving the "zero- mode constraints" (i.e., constraint equations for the longitudinal zero modes) nonperturbatively 
The zero-mode constraint appears in the discretized light-cone quantization (DLCQ) approach Q, where we set 
periodic boundary conditions for scalars in the longitudinal direction with finite extension. Of course the NJL model 
has no scalar fields as fundamental degrees of freedom, but we overcame the situation by considering the chiral 
Yukawa model. This model shows DxSB in large N limit (TV is the number of fcrmions) and goes to the NJL 
model in infinitely heavy mass limit of scalar and pseudoscalar bosons. We showed that the zero-mode constraint 
of the scalar field correctly produces a gap equation for a chiral condensate and calculated masses of the scalar and 
pseudoscalar bosons from poles of their propagators. Therefore, in Ref. we succeeded in describing indirectly the 
chiral symmetry breaking of the NJL model on the LF. 
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Since the very essence of the previous analysis was the existence of scalar fields, one may ask a question: How can 
one formulate D^SB without scalars? In order to answer the question, we treat the NJL model without introducing 
auxiliary field. An important key was already shown in Ref. 0. It was found that the "fermionic constraint" plays 
the same role as that of the zero-mode constraint. Splitting the fcrmion field as ^' = ip+ + ip± = A±'i> by 
using projectors A± = 7^7^/2, we easily find that the "bad" component ip- is a dependent variable and subject 
to a constraint equation called the "fermionic constraint." In the LF NJL model, the fermionic constraint is very 
complicated and it is difficult to solve it as an operator equation. However, we will see that to solve this equation is a 
crucial step for describing the broken phase and will find a close parallel between the fermionic constraint for DxSB 
and the zero-mode constraint for spontaneous symmetry breaking in scalar models. Although such special importance 
of the fermionic constraint might be restricted only to the LF NJL model and therefore most of the analysis might 
be model dependent, but what we are eventually interested in is the physics consequences of the chiral symmetry 
breaking. And of course we cannot reach the chiral symmetry breaking in QCD unless we understand the simpler 
and typical example of this phenomenon. Therefore the importance of our analysis is evident. 

Let us comment on other attempts for the NJL model on the LF. First of all, Heinzl et al. |^ treated the model 
within the mean-field approximation and insisted delicacy of the infrared cutoff to obtain a chiral condensate. The 
meaning and necessity of such cutoff scheme was clarified in Ref. [^. As mentioned above, an observation that a gap 
equation for a chiral condensate emerges from the fermionic constraint was first pointed out by one of the authors [[7| . 
The light-cone (LC) wavefunction of a pionic state was calculated through the LC projection of the Bethe-Salpeter 
amplitude which was derived in the equal-time quantization [|l0|. Bentz et al. introduced the auxiliary fields to 
fermion bilinears and solved the constraint equations for them by 1/iV expansion jl^. They obtained "effective" 
Lagrangian for the broken phase and discussed the structure function of the pionic state. With all these studies, 
however, there still remains many unknowns concerning basic problems. Especially, we still do not understand well 
the LF chiral transformation itself. To what extent is it different from the usual chiral symmetry? How is the chiral 
symmetry breaking realized on the LF? These fundamental problems will be resolved in the present paper. 

The paper is organized as follows. The rest of this section is devoted to introduction of the NJL model and our 
notation. In Sec. II, we discuss the complexity of the fermionic constraint in great detail. We explicitly solve the 
fermionic constraint in classical treatment and investigate properties of LF chiral transformation. In Section HI, we 
solve the fermionic constraint in quantum theory by the 1/N expansion. Here we introduce the boson expansion 
method in order to solve the bilocal fermionic constraint with systematic 1 /N expansion. We see emergence of the 
gap equation for the chiral condensate from the fermionic constraint. We obtain the Hamiltonian with respect to the 
(bilocal) bosons which is introduced by the boson expansion method. In Section IV, some physics consequences of 
the chiral symmetry breaking are discussed. First of all, we see how the chiral symmetry breaking is realized in the 
LF formalism. We discuss unusual chiral transformation of fields and nonconservation of the light-front chiral charge. 
Secondly, we construct the bound-state equation for mesonic states and solve it for scalar and pseudoscalar mesons. 
Thirdly, we derive the partially conserving axial current (PCAC) relation. Summary and conclusion are given in the 
last section. Miscellaneous topics with detailed calculation are presented in Appendices. 

Before ending this section, let us fix our model and notation. Since the primary purpose of our paper is to study 
basic properties of the LF chiral symmetry, we consider only one flavor case for simplicity. Thus the model we discuss 
is 

0^ - 

£ = *(z^-mo)*+ y [(**)2 + (*175^')2] . (1.1) 

Here ^'^(a;) (a — 1,...,N) is a four component spinor with "color" internal symmetry U{N), which has been 
introduced so that we can use the expansion as a nonperturbative technique. We always work with a nonzero 
bare mass ttiq 7^ 0. The primary reason is that the Hamiltonian with a massless fermion is plagued with a troublesome 
situation in 1-f 1 dimensions: As we will see, if we set mo = from the beginning, the canonical LF Hamiltonian P~ 
of the Gross-Neveu model vanishes at all. Even in 3-1-1 dimensions, we will see that absence of the bare mass term 
causes an inconsistency of the results. The secondary reason is to avoid massless particles which can move in parallel 
with =constant surface. The difficulty of describing massless particles is intimately connected with the fact that 
on the LF, the (massless) Nambu-Goldstone boson becomes physically meaningful only when wc first include explicit 
breaking term and then take the vanishing limit of it |p^ . The same situation was observed in the chiral Yukawa 
model |]. 

In practical calculation, it is convenient to introduce the two-component representation for the gamma matrices so 
that the projectors A± are expressed as 

A.^(; I). A-^(S I). (1.2) 
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Then, the projected fermions have only upper or lower components: 



V^+=A+vI/^2-^(^^j, V-=A_M'^2-^ (^Jj, (1.3) 

where we defined two-component spinors ip and x- Among various representations which satisfy Eq. we choose 

a representation having a similar structure to the chiral representation in 1+1 dimensions, i.e., = a^, 7^ = 
— ifj^, 75 — 7°7^ — (T^. Exphcitly, they are 

for z = 1 , 2 and 



^-li ; ■ r'=(l „ . 7'= „ ,;.]. (L4) 



Results of the chiral Gross-Neveu model in 1+1 dimensions can be easily obtained if we make a replacement for the 
Pauli matrices ct^ ^ 1 and cr' — > 0, and regard tp and x one component spinors. The explicit form of the Lagrangian 
in the two-component representation is given in Appendix A. 

In the previous work [Q, we made the longitudinal direction finite in order to carefully treat the longitudinal zero 
modes of the scalar fields. However, in the present analysis, we work in an infinite longitudinal space. There is no 
need of introducing finite . When we take the inverse of i9_, we need to specify the boundary conditions. We here 
follow the conventional antiperiodic boundary condition^ which is standard for free fermions: 

i>a{x^ ^ ~CO,x]_) ^ -Ipaix^ ^ O0,x"j_) , (1.6) 
Xa{x^ ^ ~00,x]_) ^ ~Xa{x^ = 00,x'j^) . (1.7) 

So our results must always have a smooth free field limit. 

II. COMPLEXITY OF THE FERMIONIC CONSTRAINT AND THE CHIRAL SYMMETRY 

It is highly complicated structure of the fermionic constraint that makes the analysis of the LP NJL model difficult. 
However, we cannot know anything about LP chiral symmetry unless we confront with this complexity. Therefore 
in this section, we investigate the fermionic constraint in great detail. Pirst of all, we classically solve the fermionic 
constraint. Using the explicit solution, we then discuss properties of the LP chiral transformation. Especially we show 
that LP chiral transformation is no longer an exact symmetry when mg 7^ 0. Pinally we consider the implication of 
the fermionic constraint in quantum theory. 

A. Fermionic constraint and its classical solution 

The fermionic constraint is immediately obtained as the Euler-Lagrange equation for x- 

td^Xa = ^ {-(J'd, + mo) V'a - Y {V'a (V'lxfc + xUb) + f^Va {^'^Xb " xWa) } , (2.1) 

where summation over color and spinor indices are implied. If we want to solve this equation as an operator equation 
in quantum theory, we need a commutation relation between x with ^ which must be given by the Dirac bracket. 
Since the anticommutator ip} is very complicated, it seems almost hopeless to find an exact quantum solution of 
it. However, in a classical theory where we treat all the variables just as Grassmann numbers, the equation becomes 
tractable and it is not difficult to solve it. Indeed, the exact solution with antiperiodic boundary condition is given 
by (see Appendix B for more details) 



'^For a scalar field, antiperiodic boundary condition in infinite longitudinal space leads to inconsistency [Q. However, fermionic 
fields are free from such troubles. 
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Xla{x) 



1 

71 



dy Gab{x ,y ,x±) 



mo-)jJib{y 



dz'ip2b{y ) 
-^o^lbiy~) 



where dz ^ di — 182 and the "Green function" Gab{x , y , xj_) is 



Gab{x',y-,x±)^G^''Hx) 



1 



G(°)(a;)=Pe^' 



7: 
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A{y-)dy- 



e{x -y ) + C 



G^'^-\y) 



•^ijab 



(2.2) 



(2.3) 

(2.4) 
(2.5) 



The integral constant G is determined so that the solution satisfies the antiperiodic boundary condition. In Eq. (^.4|), 
P stands for the path-ordered product. Note that we easily derive the solution for the chiral Gross-Neveu model by 
extracting the 1-1 component of A and neglecting dz- The result is equivalent to the solution of the Thirring model 
obtained by Domokos |1^. And also, if we take the free fermion limit 0, we of course recover the free solution 

due to G^°\x) 1 and G{x- ,y-) t{x- - y-)/2i. 



B. Chiral symmetry on the Hght front 



Since the bad component x is a constrained variable in the LF formalism, we impose the chiral transformation only 



on the good component ?/;+ 



■0+ or in the two-component representation [see Eq. (1.5)] 



(2.6) 



Now we have completely solved the fermionic constraint for x, we can explicitly demonstrate its transformation 
property under the LF chiral transformation. However, before discussing the NJL model, it will be instructive to 
remind you of the LF chiral symmetry in the free massive fermion. 



As we mentioned before, the massive free fermion is chiral invariant under the transformation (2.6). Let us see 
this fact directly in the Lagrangian even though it is a little lengthy. It is convenient to separate the solution of the 
fermionic constraint x — {V^id-)^^ (^—a^di + toq) tp into mass-independent and dependent parts x = X^^^ + x*-™-* as 



\/2 id- 



Note that there is a relation between x*-"-* and x^"^''- 

a'd.x^""^ + moX*°^ = 



(2.7) 



As a result of the LF chiral transformation 



we find 



X 



(m) 



/(O) 



(2.8) 
(2.9) 



X^(^FC, where weom = id+ijj T^icr^d, 



V2^ 

is the fermionic constraint. The 



The free fermion Lagrangian is compactly expressed as £froc = 

i^o)x = is the equation of motion for tjj and wpc = id-X — + "t-o)V' 

second term is zero and is invariant under the LF chiral transformation. Now substituting x — X^'^'' + x''™^ vato the 
Lagrangian, the first term is decomposed into apparently invariant and (seemingly) non-invariant terms 



m^X 



(m) 



^(.'ftx'"" 



m^x 



(0) 



The first term consists of the mo-independent term and quadratically dependent term 0{m'^), while the seco nd t erm 
linearly depends on mp. The O(mo) term changes under the chiral transformation, but due to the relation (2.7), it 
eventually vanishes and therefore the Lagrangian is invariant even if there is a mass term. As a result, we have a 
conserved Noether current 101 



4 



^M.?5Froc — I 

which of course reduces to the usual current in the massless hmit. 

Now let us consider the NJL model. Decomposition of x is straightforward: 



=-^y_rf2/-G.,(.-,,-,.,)(^— f ), (2.10) 



Since the matrix A, an d thu s Ggb jx, y) is invariant under the transformation (2.6), it is ea sy t o find that and 



x'-'"-' transform as Eqs. (2^) and (|2.9| ). Therefore, if mo = 0, the LF chiral transformation (|2.6| ) is equivalent to the 



usual chiral transformation. The chiral current and the chiral charge are given by 

35 = *7^75* , (2.12) 
dx^£x^j+{x) = / dx~(fx±ip^a^i; . (2.13) 



)LF ^ / j^-j2„ 
^ — oo 



How about the massive case? As we explicitly showed above, the mass term does not prevent chiral symmetry in the 
free fermion case. We must bare in mind such possibility even in the NJL model. Thus it is worth while to check 
whether the massive NJL model is invariant under the LF chiral transformation. To see this, it is convenient to treat 
the Hermite Lagrangian 

>CHcrmitc = 5+ V - ^ [(^V'S^ + ^.X^V) + "iQ (V^X + xV)] • (2-14) 

Note that this is equivalent to the free Lagrangian except that x is a solution of Eq. (^.l|). Now the apparently 
non-invariant term is a term linearly depending on toq: 

- (fT^Sa^'"^ + mox^"^) + H.c. . (2.15) 

In the massive free fermion case, we had the same term but it eventually vanished due to Eq. (|]^). However, in 
the NJL model, it is evident from Eqs. (2.10) and ( |2.11 ) such relation does not hold because G depends on xj_. 



Therefore we have verified that the massive NJL model is not invariant under the LF chiral transformation. If and 
only if mo — 0, the LF chiral transformation is the symmetry of the NJL model and equivalent to the usual chiral 
transformation. This is of course not a surprising result but must be checked explicitly. Anyway, we do not stick to 
this problem anymore. 



Irrespective of whether we have a mass term or not, we always use the definition for the chiral current Eq. (2.12) 
which was derived for the massless fermion. In the massless case, it is of course a conserved current 9pj5 = 0, while 
in the massive case, a usual relation holds 

d^j^ = 2mo'l'i75* , (2.16) 



which is derived by using the Euler-Lagrange equation for the massive fermion. Equation (2.16) is used when we 
discuss nonconservation of the chiral charge and the PCAC relation in Sec. IV. 



C. Implication of the fermionic constraint 

So far we treated the fermionic constraint in classical theory and obtained the exact solution (|]^). However, this 
solution does not give a nonzero condensate and the resulting Hamiltonian does not describe the broken phase. The 
situation is very similar to the previous analysis of the chiral Yukawa model Q| . The chiral Yukawa model in the 
DLCQ approach has three constraint equations. We solved them in classical theory but we could not find any way 
to describe the broken phase with the classical solutions. What we finally found is that it is very important to treat 
the constraint equations, especially the zero-mode constraint for a scalar field, nonperturbatively in quantum theory. 
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This fact is true of our present case. To obtain a nonzero condensate, we must treat the fermionic constraint as an 
operator equation and solve it with some nonperturbative method. 

To strengthen this, let us briefly overview the procedure in the previous work [Q. In the chiral Yukawa model, 
we have three dependent variables: Two of them are the longitudinal zero-modes of scalar and pseudoscalar fields 
o'oix^) — (1/2L) J^j^dx~a{x), 7ro(a;_L) = (1/2L) J^j^dx~Tr(x) where L is an extension of the longitudinal direction 
x^ e [~L,L], and the rest is the bad component of a fermion So there are three constraints: 



N V2 



dx 
'2L 



«75 



at 



-1 



(2.17) 



2i9_V'° = {il\_d, + mo + (T - ^7175) 7+^+ , (2.18) 

where A — g^N in the present notation and /i is a dimensionless parameter which controls the scalar and pseudoscalar 
masses. In the infinitely heavy mass limit, fj, 00, we re cover the NJL model. The procedure of Ref. [0 is as follows: 



First, we formally solved the fermionic constraint (2.18) and substitute the solution into the zero-mode constraints 
(2.17). Second, we solved the zero- mode constraints by 1/A^ expansion with a fixed operator ordering and found that 
the leading order of the scalar zero-mode constraint can be identified with the gap equation. Selecting a nontrivial 
solution of the gap equation, we again substitute it back to the fermionic constraint. Then we obtain the final 
expression for the bad component ■(/'- in terms of independent variables. Thus we solved three coupled equations 
step by step. On the other hand, we have only one constraint equation. The procedure in the chiral Yukawa model 
suggests that we will have to do almost the same procedure at once when we solve the fermionic constraint. Note 
that just the same as in the chiral Yukawa model, a perturbative solution cannot reach the broken phase even in 



quantum theory. Therefore, we naturally expect that solving the fermionic constraint (2.1) in quantum theory using 
some nonperturbative method is necessary for describing the chiral symmetry breaking . 



III. SOLVING THE FERMIONIC CONSTRAINT BY 1/N EXPANSION IN QUANTUM THEORY 



As we discussed above, it is important to solve the fermionic constraint (2.1) in quantum theory by some nonper- 
turbative method. Here we solve it with a fixed operator ordering by using the 1 /N expansion. For systematic 1 /N 
expansion, we introduce the bilocal formulation. We rewrite the fermionic constraint in terms of bilocal fields and 
expand it following the Holstein-Primakoff type expansion of the boson expansion method. We always work with 
fixed x^ . 



A. Quantization and the operator ordering 



To solve the constraint in quantum theory, we must first perform the Dirac quantization for constraint systemsQ. 
After tedious but straightforward calculation of the Dirac brackets, we find a familiar relation for the good component 
i>a {a = 1,2) 



lew, 4^iy)} ^ ^ = So^pSabSix- - y-)5'^'Hx^ - y^) 



and so on. Wc introduce the simplest mode expansion at x"^ — as in Ref. |14[: 



d^ki 

~2^ 



dk+ 
V2^ 



(3.1) 



(3.2) 



where kx = —k~^x + fc^a;*^. The vacuum is defined by the annihilation operators as 



^ In Ref. the authors solved the constraint equations for auxiliary fields before canonical quantization was specified and 
gave a c-number to the scalar auxiliary field in leading order of 1/A'^. Nevertheless, the condensation in the NJL model is a 
quantum phenomenon and thus this procedure is not justified. 
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bl{k)\0)=dl{k)\0)=0. (3.3) 

When we deal with the quantum fcrmionic constraint, we have to specify the operator ordering. In many pubhcations 
discussing the zero-mode constraints, people often choose the Weyl ordering with respect to both constrained and un- 
constrained variables. However, in a previous paper we discussed that the ideal situation was to find a "consistent" 
operator ordering. For example, let us consider an anticommutator {x, tp} in the NJL model. It can be evaluated in 
two different ways: (i) by using the solution Xsoi = xii') of the fermionic constraint and the standard quantization 



rule (3.1), and (ii) by calculating the Dirac bracket for {x, ''/'}• For the case (i), we select a specific operator ordering 
for the fermionic constraint, and the result depends on the ordering. For the case (ii), we must also determine the 
ordering in the r.h.s. of the Dirac bracket {Xj'^Id = • ■ ■• These two results must be equivalent to each other. We 
have two ambiguity of the operator ordering: those of the constraint equation in (i) and the right hand side of the 
Dirac bracket in (ii). "Consistent operator ordering" should be imposed so that these two quantities be identical. 
In other words, we determine the operator ordering of the r.h.s. in the Dirac brackets so that it coincides with the 
direct evaluation. In the chiral Yukawa model, we could not check that the ordering we adopted was consistent or no t. 
Again in the NJL model, this is a very difficult task and we choose a specific operator ordering defined by Eq. (^.]|). 
However, the chiral Gross-Neveu model in dimensions allows us to check the consistency of this operator ordering. 
This is briefly shown in Appendix C. 



B. Boson expansion method as 1/N expansion of bilocal operators 



How can we solve the "operator equation" Eq. (2.1) by the 1/A^ expansion? It is generally difficult to count the 
order 0(iV") of an operator instead of its matrix element. What is worse, it is physically hard to justify the 1/A^ 
expansion of the fermionic field itself. However, as was discussed in Ref. [|l^, there is a powerful method to this 
problem. We can perform a systematic 1/A^ expansion of operators if we introduce the bilocal operators and use 
the boson expansion method. The boson expansion method is one of the traditional techniques in nonrelativistic 
many-body problems |]l6| . Originally this was invented for describing bosonic excitations in non-bosonic systems such 
as collective excitation in nuclei or spin s yste ms. 

Let us rewrite the fermionic constraint (2J) in terms of bilocal operators. We introduce the following "color" singlet 
bilocal operators at equal light-front time 



N 



, N 
. N 



(3.4) 
(3.5) 
(3.6) 



We define the Fourier transformation of them as 



[°° d'^x f°° d^y 
M^,ip, q) = M^p{x, y 



^ ipx iqy 



and so on. Note that this definitio n al lows the longitudinal momenta to take negative values. Using these bilocal 
operators, the fermionic constraint (2.1) is cquivalcntly rewritten as 

i-^Taf}{x, y) = ^ {up^Ma-^ix, y) - a\f^Mjaiy, x)) + ma (^Mapix, y) - Mf3a{y, xfj | 

^Ma^{x,y){5^pT{y,y) + i(j^p U{y,y)^ 

-(^5p-^T{y,y)-taj^U{y,y)^M^a{y,x)^ , (3.7) 



' 2 



and 
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'-^Uap{x,y) i ^-d^^ {al^Matix.y) + a\pM^a{y,x)^ + mo{Map{x,y) + Mpa{y,x)^^ 

^Maj{x,y)(SypT{y,y) + ia^p U{y,y)^ 

+ {5p^T{y,y) - iafj^ U{y,y)'jMya{y,x)^ , 



2 



(3.8) 



where we have introduced quantities T{x,y) = Taa{x,y) and U{x,y) = {<J^)ap l^ap{x,y) so that ^'^{x) — T{x,x) 
and 4'«75^'(a;) = U{x,x). In actual calculation, it is more convenient to treat equations for the operators without 
spinor structure T{x,y) and U{x,y) because they form closed equations (see Appendix D). Once we solve them, we 
immediately obtain Tap{x, y) and Uap{x, y) from the above equations. 

For systematic 1/N expansion of the bilocal fermionic constraints, one must know how to expand Map{p, q)- It is 
the boson expansion method, especially, the Holstcin-Primakoff type expansion for large N theories, that enables us 
to expand Mapip^ q) as operator quantities: 



mMp,q)-nY,[-7^) &p,q)- 



(3.9) 



According to the Holstein-Primakoff expansion [Eqs. (D5)^(D8)], the first three terms are written in terms of bilocal 
bosonic variable B{p, q) as 



f^''^iip,q)=Sc.(3S^'Hp + q)9{p+)0{-q+) , 

^^(^,9) = Bp^{q,p)e{p+)e{q+) + Blj{-p,-q)e{-p+)e{^q+) , 
mSIp,?) = / [dk]Y,Bi^{-p,k)Bp.,iq,k)9{-p+)0{q+) 

- j [dk]J2Blp{k,^q)B,^{k,p)0{p+)e{-q+) . 



(3.10) 
(3.11) 



(3.12) 



where 



[dq] 



dq^ 



Any commutator between Map{p,q)''s [such as Eq. (|D3|)] is correctly reproduced if one uses the following bosonic 
commutators 



Bai3{Pi,P2), B\g{q^,q2) ^ So^^SpsS^^^p^ - q^)S^^'^ (p^ - q2) , 
BapiPi:P2): B^s{qi,q2) =0 (p+,g+>0). 
Note also that the state annihilated by B{p, q) is identified with the original Fock vacuum: 

B{p,q) |0) = 0. 



(3.13) 
(3.14) 

(3.15) 



More detailed discussions about the boson expansion method applied to LF field theories are found in Ref. ||l5| and 
Appendix D of the present paper. 



C. Solution to the bilocal fermionic constraint 



We are ready to solve the bilocal fermionic constraint using the 1 A/V expansion. As we commented before, it is 
convenient to solve the equations for T{p,q) and U{p,q) (see Eqs. (Dl) and ( ^^ ) in Appendix D for their explicit 
forms). Once we know T(p, q) and U{p, q), then it is straightforward to obtain Tap{p, q) and Uap{p, q)- 

Expanding T(p, q) and U{p, q) as 
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r(p,g) = A^^(-=j t^"Hp,q), (3.16) 

n— ^ 

oo . s n 

W(p,g) = 7V^(-=j (3.17) 

and inserting them into the fermionic constraints, we find for the lowest order 0{N) 

«(°)(p,g)j~l 'j-5<^^i_(2^U(°)(fc,p + q-fc)j' 

where ijg = 5^-/V. Since there are no operators in these equations, i^*'-' and u^^^ are c-numbers. Nonzero solutions give 
leading order contribution to {^'^) and (^275^), 

(2^ 

(0| ^*75* m^Nj -^Mp) + ■■■ , (3.20) 

where t^'^''{p,q) — ft{p)S^^\p + (l) and u^'^^p^q) = fu{p)5'^^'' {p + q). As the above equation ( |3.1j ) with toq = is 
invariant under the chiral rotation, we can always take u^*''(p, q) — 0. For the massive case, we also take it(°^(p, q) = 
and t^"^ (p, q) ^ 0. Now let us introduce a quantity M, which corresponds to the dynamical mass of fermion, 

A/ = mo-5^(§*) . (3.21) 

Then, to obtain t^^\p,q) is equivalent to determining M, viz. 

tio}(p,q) = -M'-^6^'\p + q). (3.22) 



(0|^*|0)=iV / j^ftip) + ••• , (3.19) 



In terms of Af , the leading order fermionic constraint (3.18) is rewritten as 



M-mo 2 f (Ppe(j)+) 

^ .9o / TTT-v?— ^ • (3-23) 



M -'"J (27r)3 p 

Physically this equation should be interpreted as a gap equation. This is clarified in the next subsection. 

Similarly higher order fermionic constraints are solved order by order. This is because the fermionic constraints for 
t'^")(p, q) and w'^"-'(p, q) are linear equations with respect to the highest order. The n = 1,2 solutions are important 
for giving a nontrivial Hamiltonian and so on. More details are discussed in Appendix D. 



D. Gap equation 



Now let us discuss the physics meaning of Eq. (3.23). As we mentioned above, this equation should be regarded as 
a gap equation for chiral condensate. In several previous studies of ours, we have seen essentially the same kind of 
equations Indeed, in Ref. 0| it was pointed out that Eq. ( 3.23D itself is the gap equation. Also in the chiral 

Yukawa model |H] , the zero-mode constraint for the scalar field reduced to the above equation and was interpreted as 
a gap equation. Since this identification is an indispensable step for our framework, let us again explain it within the 
NJL model. 

First of all, consider a naive massless limit toq ^ of Eq. ( 3.23 ): 



M 1 



9o 



(fp e(p+) 
(27r)3 p+ 



Thus we find two possibilities: the first is M 



9l 



and the second is 
(Pp e(p+) 



(27r)3 



= 



(3.24) 
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but M is arbitrary. Of course the first case is not interesting because it corresponds to the symmetric phase. On the 
other hand, the second case with nonzero M does not immediately mean the existence of broken phase. Since the 



equation (3.24) is independent of M as it is, the dynamical mass M is left undetermined, which is not a ph ysically 



acceptable situation. However, this obse rvati on is not correct because the divergent integral in Eq. ( 3.24 ) is not 



regularized. Indeed, we can identify Eq. (3.24) with the gap equation only after we carefully treat the infrared (IR) 
divergence. 

To see this, let us put an IR cutoff. First consider the same cutoff schemes as in the equal-time formulation, such 
as the covariant four-momentum cutoff. We can easily translate it into a cutoff on the light-cone momentum fc+ and 
and obtain the correct gap equation. Indeed, in Ref. a noncovariant (rotationally invariant) three-momentum 
cutoff was performed to obtain the known result. But such a cutoff is artificial as a light-front theory, and we here 
adopt another cutoff scheme, the parity invariant cutoff. Usually, it is natural and desirable to choose a cutoff so as 
to preserve symmetry of a system as much as possible. For the LF coordinates and x\, it would be natural to 
consider parity transformation (a:"*" <-> a;^ — > —x\) and two-dimensional rotation in the transverse plane. In the 
usual canonical formulation where a:"*" is treated separately, the parity invariance is not manifest. However, we find it 
useful for obtaining the gap equation. In momentum space, the parity transformation is exchange of A:+ and fc^ and 
replacement k^^^ — > —k\. Therefore the parity invariant cutoff is given by fc* < A and k\ < A'^. Using the dispersion 
relation 2fc+fc^ ^ ^± ~ -^^^j fi^d that the parity invariant regularization inevitably relates the ultraviolet (UV) 
and IR cutoffs: 

+ <k+ <A. (3.25) 
2A ^ ^ 

This also implies the planar rotational invariance k^ < 2A^ — = A'^. What is important here is the use of 
constituent mass M in the dispersion relations. Phy sically it corresponds to imposing self-consistency conditions. 



Since the IR cutoff includes M, the r.h.s. of Eq. ( 3.23 ) has nontrivial dependence on M: 



M Att^ \ A2 V M2 J 

This is the gap equation and is equivalent to that of the previous result in the chiral Yu kawa model [Q . It has a 
nonzero solution M ^ even in the mg limit. The somewhat unfamiliar equation ( 3.26| ) of the NJL model 



exhibits the same property as the standard gap equations of the equal-time quantization. For example, when mo = 0, 
there is a critical coupling g'^^. = 2tt'^ /A"^, above which M 7^ 0. 

The essential and inevitable step to obtain the gap equation is the inclusion of mass information as the regularization 
rather than the fact that the UV and IR cutoffs are related to each other. If we regulate the divergent integral 
without mass information ( e.g., introducing the UV and IR cutoffs independently), we cannot reproduce the gap 
equation. The loss of mass information is closely related to the fundamental problem of the LF formalism [ p7[ , and 
the parity invariant regularization can be considered as one of the prescriptions for it. Reference |]l7| discussed within 
scalar theory that the light-front quantization gives a mass-independent two-point function (at equal LF time), which 
contradicts the result from general arguments concerning the spectral representation. We have been encountered with 



the same problem in Eq. (3.23) because the integral is regarded as a naive estimation of (\E'^')/M by using fermion 



with mass M. And also the origin of mass-independent result can be traced back to the mode expansion (3.2). Even if 
we include the wave function for free fermion field, we do not have any mass-dependence on the mode expansion [ p^ . 

Let us give a brief comment on the chiral Gross-Neveu model. Of course the important difference of the H-1 
dimensional case is the renormalizability. Ignoring the transverse directions in the above calculation, we easily find 
the gap equation (M — mo)/M = go/(27r) \n{2A'^ /M^) where the parity-invariant cutoff AP/2A < A:+ < A was used. 
Though it explicitly depends on the cutoff A and is divergent as A — > cx), we can remove the divergence by coupling 
constant renormalization 0|. 



E. Hamiltonian 



Having the solution to the bilocal fermionic constraint, we can rewrite the fermion bilinear operators in terms of 
the bilocal bosons. Of special importance is the (Hermitian) Hamiltonian, which is easily expressed by 7^^(p, q) and 
Uaj3{p,q) as follows 

H = P- = ^ I d^x [(V'V^S.x + da^<j'^) + nio (V^X + xV)] 
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mo 



Jd^p d\T{p,q) 6^^\p + q) 



(3.27) 



Apparently this Hcrmitian version of the Hamiltonian seems equivalent to the free Hamiltonian, but the information 
of interaction enters through the bad component x- We find the 1/iV expansion of the Hamiltonian 



H = N 



OO 

E 

n=0 



(3.28) 



by substituting the solutions of the fermionic constraints into the Hamiltonian. The zeroth order contribution is just 
a divergent constant and we discard it. The first order is strictly zero. Nontrivial contribution comes from the order 



J[dp][dq] 



bIAP' Q)Bai3{p, q) 



2p+ ' 2q+ J 
dp][dq][dk][dl]S^^'>{p + q-k- l)a{p+ 



q 



(27r)3 

X [5^^f(p,g;fc,Z)-I]"f(P,g;fc,0] Blp{p,q)B^s{k,l) 
c number , 



where a{p'^ +9^) is defined by Eq. ( |D13| ) and "kernels" of the interaction terms are 



J2 ;f (p, g; k, I) ^ \s{-p) - si-q)] \s{k) - s{i) 

— L J aP - 



V{-p)-V{q) V{~k)~V{V} 



a/3 



(57 
^7 



with 



Sa0{p) 



ip^a"^ 



M 



2p^ 



Vo.f}{p) 



Q/3 



ip'a' - M 3 
2p+ 



(3.29) 

(3.30) 
(3.31) 

(3.32) 



Q/3 



As is evident from the explicit forms of the kernels ( p. 30 ) and (3.31), they originate from the scalar interaction (4'^')-^ 
and the pseudoscalar one (^'175 4')^, respectively. If we substitute a nontrivial (trivial) solution of the gap equation 
(3.26) into the above Hamiltonian, then it governs dynamics of the broken (symmetric) phase. The first term of /i^^^ 
corresponds to a free part with mass M and the second term to an interaction part. In the broken phase, M is the 
dynamical mass and the Hamiltonian suggests a constituent picture. 

As we mentioned before, the Hermite Hamil tonia n of the chiral Gross-Neveu model has only an mo-dependent term. 
Neglecting the transverse coordinates in Eq. (3.27), we have 



GN 



Too 
2^/2 



dx (V'^X + xV) 



Furthermore, the classical solution for the bad spinor component x is proportional to mo. Therefore naive mo 
limit gives just a zero Hamiltonian. However, if we solve the gap equation and substitute the nontrivial solution into the 
Hamiltonian, the resulting Hamiltonian turns out to be prop ortional to and survives even in the chiral limit. This 
is easily seen from the Hamiltonian of the NJL model ( 3.29 ). The (constituent) mass term in Eq. (j3.29| ) comes from 
the bare mass term, whose factor mo cancels with a factor M^/mo in the second order solution / d'^p7^\p, —p). Of 
course this is not reached if we set mo = from the beginning. Therefore inclusion of the bare mass term is necessary 
to obtain a correct (constituent) mass term of the Hamiltonian. 



IV. PHYSICS IN THE BROKEN PHASE 

By solving the fermionic constraint, we acquired necessary ingredients for discussing physics consequences of the 
chiral symmetry breaking. Basic quantities such as ^'175^', and the null-plane chiral charge (2.13) are expressed 
in terms of the bilocal bosons Baj3{p, q) and B^^p{p, q) as 
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N 

^^(x) = T(x,x) = -^(ma - M) 
9o 



(27r)3 



(27r)3 
u^^\p,q) e*(P+g)^ 



t^^\p,q) e^(P+9)^ + 0(7V") , 
0(iV°) , 



(4.1) 
(4.2) 
(4.3) 



where t'^^'^ij)^ q) and u'-^-'(p, q) are given in Appendix D. Now that thes e are given as fu nctions of the bilocal bosons 
at the operator level, all the calculation is done with the commutators (3.13) and (3.14). 



A. Chiral transformation and nonconservation of chiral charge 

Why could we obtain a nonzero fermion condensate? To understand this, let us rewrite the fermionic constraint 



as 

^2 



id^Xa = {-cr'di + mo) tpa - {tpaT{x, x) + ia^4'all{x, x)) 



and substitute Eqs. (4.1) and ( |4.2| ) into T{x,x) and U{x,x), respectively^. Then the leading order equation turns 



out to be equivalent to the constraint equation for a free fermion with mass M, 

V2 

Also at the same order, the equation of motion for the good component ip says that the fermion acquires a mass M. 
This means that the operator structure of the bad spinor x changes depending on which solutions of the gap equation 
( |3.23 ) is selected. For massive fermion, the fermion condensate (^f^*) is no longer zero even if the vacuum is trivial. 



One can find an analogy between the chiral Yukawa model and the NJL model because in the chiral Yukawa model, 
operator structure of the longitudinal zero modes and subsequently of the bad spinor component changes depending 
on the phases. 

One thing to be noted is the peculiarity of the mode expansion (|3.2| ). It is evident that the mode expansion 
has no mass dependence in it. This caused the problem of identifying the lowest fermionic constraint with the gap 
equation. We had to supply mass information properly when we regularize the IR divergence. On the other hand, 
such independence of mass, in turn, implies that our mode expansion allows fermions with any value of mass. In other 
words, the LF vacuum does not distinguish mass of the fermion. Therefore we can regard the vacuum for massless 
fermion as that for massive one. Mass of the fields is determined by the Hamiltonian. This is the reason why we can 
live with the trivial vacuum while having a nonzero fermion condensate. This fact is not a limited phenomenon for 
our specific mode expansion but a common one for light-front field theories. Indeed, even if we expand a fermion field 
with free spinor wave functions, u(p) and v{p), we have no mass dependence | |To| . 

The fact that the operator structure changes depending on the phases, also resolves a seeming contradiction between 
the triviality of the null-plane chiral charge and the nonzero chiral condensate (0| |0) 7^ 0. In general, it is known 
that a null-plane charge always annihilates the vacuum irrespective of the presence of symmetry. This can be checked 



explicitly by the expression (4.3), viz 



|0) = . (4.4) 

However, the triviality of Q]^^ in the presence of the chiral condensate immediately leads to a contradiction if an 
equation [(55^,^*75^] = — 2i4"I' could hold in the broken phase. In the previous analysis of the chiral Yukawa 
model [Q , we were faced with exactly the same problem and resolved it by recognizing that in the broken phase the 
chiral transformation of dependent variables are different from the usual one simply because their operator structure 



^In the leading order, this procedure corresponds to the mean-field approximation done by Heinzl et al. They solved the 
fermionic constraint by simply linearizing the interaction parts as —g^/^/2'>jja{T{x,x)). By evaluating the vacuum expectation 
value {T{x, x)) self-consistently with the dynamical fermion mass M = mo — g^{T{x, x)), they obtained the gap equation. If 
one uses the parity-invariant cutoff, the result coincides with ours. 
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changes. This is of course true of the NJL model. First of all, as we saw above, if we select the nontrivial solution 
of the gap equation, the fermion is no longer a massless fermion even in the chiral limit. Secondly, we can explicitly 
show that the usual transformation law [Ql"^, ^175^] = — 2i^'^' holds only in the symmetric phase {M = 0). In the 
broken phase a simple calculation (up to 0{N^/'^)) leads to 



N 



= -2i-9^{x) + 2i — (too - M) + 2iVNM 
9o 



(1)/ ^ 



2lb^ 
■90 ^+ 



dpdq 



iip+q)x 



) / _|_ , A--. [ d k fiao 



(27r)3 ^ 
d^k ^al{k,p + q- k) 



k+ 



-0{N°) 



(4.5) 



Even if we take the chiral limit toq — > 0, the extra term proportional to M survives nonzero. This also implies that if 
we put M = 0, the usual relation holds. The unusual chiral transformation, however, is consistent with the triviality 
of Q^^ because (0| [Q^'',§i75*] |0) = 0. 

A similar situation occurs for the Hamiltonian. Nonconservation of the null-plane chiral charge has been pointed 
out by several people as a characteristic feature of the chiral symmetry breaking on the LF . They discussed it 

under the assumption of the PCAC relation, but we can check it explicitly by using the broken Hamiltonian. After 
lengthy calculation, we find the commutator [Qs^, H] is really nonzero and again proportional to the dynamical mass 
M: 

[QLF,i/] =M^^ J [dp][dq][dk][dl]S^'Hk + l-p-q)aip++q+)Bl^ip,q)B^s{k,l) 



'Pl_9l 
q+ 



'a(3 



k+ 



1 

1+ 



^3 



k+ 1+ 



1 

5s^ 



1 



k+ 1 + 



""Si 



(4.6) 



Therefore, the LF chiral charge is not conserved even in the chiral limit. In our framework it would be more 
understandable to mention that the Ha milto nian is not invariant under the LF chiral transformation in the broken 
phase. The broken phase Hamiltonian ( |3.29 ) has three terms: Af-independent, linearly dependent and quadratically 
dependent terms. The quadratically dependent term, as well as the M-independent one, does not break the LF chiral 
symmetry. It is the term proportional to the dynamical fermion mass M which breaks the LF chiral symmetry. And 
also, since Eq. (L6) is proportional to g^, the symmetry breaking term purel y com es from the interaction^. 

This result should be consistent with the current divergence relation Eq. ( 2.16| ). Integrating it over the space, we 
have 



a 



2mo / dx d x±^ij5^ . 



(4.7) 



Therefore if the LF chiral charge is not conserved in the chiral limit, the r.h.s. must show a singular behavior 



dx d xj^^'i75VE' oc 



1 

Too 



(4.8) 



This can be verified directly by using the solution of the fermionic constraint. Indeed we find that J dx d^x_\_'9ij5'^ — 
J dp u^^\p, —p) is proportional to M/ruQ and gives exactly the same result as Eq. (L6). Importance of such singular 
behavior for making the Nambu-Goldstone boson meaningful was stressed by Tsujimaru et al. in scalar theories [ p^ . 
Assuming the PCAC relation, they showed that the zero mode of the Nambu-Goldstone boson has a singularity 
of 

™NG where tong is an explicit symmetry-breaking mass. Our result ( p^ [) is consistent with theirs because the 
operator ^i^f^'^ is directly related to the Na mbu- Goldstone boson. Later, we will prove that the PCAC relation is 
derived from the current divergence relation ( 2.16| ). 



*For a massive free fermion, we have [Q\^ ,H] — 0. 
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B. LF bound-state equations for mesons and their solutions 



1. Single bosonic state as a fermion-antifermion state 

In our formulation with the boson expansion method, any bosonic excited state is described by the Fock states of 
the bilocal bosons constructed on the trivial vacuum: 

n<ft(p.,g.) |0) . (4.9) 



Since the Hamiltonian (3.29) is quadratic with respect to the bilocal bosons, the first excited state is given by a 
single bosonic state. In fermionic degrees of freedom, the one boson state corresponds to the leading contribution (of 
expansion) of a fermion-antifermion state. To see this, let us write a mesonic state only with a "color" singlet 
fermion-antifermion Fock component: 

1 P°° 

I meson;P+,Pj_) = ^ / dk+ d?k^ (k)W^ {k)(f^ {P - k) \{)) , (4.10) 

VN Jo J-oo 

where the LC wavefunction $"'^(A;) is normalized so as to satisfy the condition 

( meson; P \ meson; Q ) = {27rf 2 P+d'-^\P - Q) , (4.11) 

or equivalently, 

•'^ a/3 



According to the Holstein-Primakoff type expansion (DS), the fermion-antifermion operator fej^d^j can be equivalently 
rewritten as 

b^J{k)dfiP - fc) = : M-^i-k, -P + k): 
= VNBl,{k,P-k) 

J[dq][dq']Bl^{q,P-k)BUk,q')B^5iq,q') + --- ■ (4.13) 
Therefore at the leading order of expansion, the mesonic state is described as a single (bilocal) boson state, 



rP^ roc 

i;P+,P±)^ dk+ d^k^^"'^{k)Blp{k,P~k)\Q)+0{N-^/^) 

Jo J-oo 



(4.14) 



Besides this, it is evident from the normalization condition ( 4.11 ), a local operator a^{P) = / d'^k^°'^{k)Bl^„{k, P — k) 



satisfies the usual bosonic commutators. 

The LC wavefunction ^"'^(jt.) ^^^^ ^he mass of a meson Mmcson is determined by solving the light-front eigen- value 
equation: 

I meson; P+, Pj_ = 0) = \ meson; P+ , Pj_ = ) , (4.15) 

where we set P_j_ — 0, for simplicity. 

2. Scalar and Pseudoscalar mesons 
In the leading order of 1 /N expansion, the Hamiltonian has only quadratic terms of the bosonic operators. Therefore, 



diagonalization of the Hamiltonian, or equivalently, solving the light-cone bound state equation (4.15) is straightfor- 
ward. First of all, if one notices the orthogonal property (P(— fe) — P(0)a/3 {^{^) ~ ^{^))i3a = ^ where k = {xP^ , k\^) 
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and Z = P — fe = ((1 — x)P+, —k]_), one can easily find the spinor structure for scalar (cr) and pseudoscalar (tt) states 
should be 



a ;P+,Pi = 0) =P^ 



dx / d^k^ Mx,kl){ikla' + il-2x)M} Bl^ikJ) |0) . 



(4.16) 
(4.17) 



These two states are orthogonal to each other. S omew hat nonstandard spinor structure of the m eson ic states is due 
to our specific choice of the mode expansion Eq. ( [3.2[ ) and the representation for 7 matrices Eq. ( |1.4| ). For example, 
if one rewrites the pseudoscalar field ^ij^'i/ in terms of the bilocal bosons, one finds the same spinor structure as 
that of Eq. (4.16). Note also that {75(71^^ + Af)}^^ — {(i/c^cr* + Af) c'^}^^ for a, /? = 1, 2 in our two-component 
representation for the 7 matrices. 

Spinor independent parts of the LC wavefunctions 4>T^^a{x,k\) are given as solutions of the following integral 
equations: 



1 



x{l — x) 

U2 , li,f2 



TO^0a-(x, k\) ^ 4 0cr(a:, fc^) ^"^^ 



(27r)3 a:(l - a:) 7o 
1 



a;(l — x) 



(2^)3 x(l - a:) Jo 



- (1 - 2y)HP 

y(i - y) 



Here the factor a — q;(P+) defined by Eq. (D13) is given as a result of the gap equation, 

-1 



mo 



d\± 



1 dx\ 
X ) 



M (27r)3 

Since these integral equations are separable ones, solutions are easily found 

gl M 1 



b-^{x, kj_) — — Cjr 



(jjaix, k]_) 



-c„ 



(27r)3 mo x{l ~ x) ~ {k^ + A'P)/ml 
gl M /ml-4AP 



(27r)3 mo 



1 



x{l -x)~ {kl + NP)lml ' 



(4.18) 
(4.19) 

(4.20) 

(4.21) 
(4.22) 



where and Co- are constants C-jr^a — Jq dx J d^kj_(j)jr^a{x, kl). 

Equations for m-^ and ma are derived from the normalization condition for the LC wavefunctions, viz. 



l = 5o 



M 




f d'^kjL 






ml 




mo 


(27r)3 


kl- 


1- Af2 


— m'^x{l — 


x) 


M 




f d^kj_ 




ml 


-4A/2 




mo 


(2^)3 


kl- 


M/2 


— mlx{l — 


x) 



(4.23) 
(4.24) 



These are exactly equivalent to the corresponding equations in the previous work on the chiral Yukawa model [Q, 
where we obtained them by calculating pole masses of the scalar and pseudoscalar bosons. If one uses the same cutoff 
scheme (extended parity-invariant cutoff) as in Ref. [Q: 



kl + A/2 kl + A/2 



< 2A2 



(4.25) 



the pseudoscalar mass for small bare mass mo is 



ml='^Z^ + 0{ml) , 



9tM 



(4.26) 



where a cutoff dependent factor 



1 /I + v/l-2A/2/A2\ V1-2A/2/A2 
■ In ' I V 



87r2 



1 - v/l-2A/2/A2 



47r2 



(4.27) 
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is related to normalization of a pseudoscalar state [see Eq. (4.35)]. Clearly vanish es in the chiral limit toq ~^ 
and the pseudoscalar state is identified with the Nambu-Goldstone boson. In Eq. ( 4.18|) , the first term corresponds to 
a kinetic energy part of the fermion and antifermion with the constituent mass M and the second term, a potential 
energy part. The masslessness of the pseudoscalar state in the chiral limit is re alized by the exact cancellation between 
the kinetic energy and the potential energy. Indeed, if we integrate Eq. (4.18) over x and k]_, we find 



dx / d kj_(j)T,{x,k]_) = 1- 



9„2 f 1 dx r <Pq± 



f"0 I 9„2 f 1 dx [■ d^q_L I I 
W + ^ffo Jo — J j2Fj^ J "'O 

(too ^ 0) . 



dy 



2/(1 - y) 



Therefore, m^r = is fulfilled in the chiral limit even though the fermion has the constituent mass. 
On the other hand, squared mass of the scalar state for small too is 

ml = 4A/2 + ©(mo) . 



(4.28) 



(4.29) 



At a first glance, the result to^ = 2M in the chiral limit seems to suggest a static picture of a fermion and an 
antifermion, but actually the mass of the scalar meson comes from a part of the potential energy. The kinetic energy 
cancels with t he re st of t he po tential energy. 



Equations ( 4.21 ) and ( 4.22 ) have the same functional form with respect to the variables x and fc^. However, the 



difference between to^ and m „ gr eatly affects the shape of the LC wavefunctions. This is most remarkable in the 
chiral limit: As too 0, Eq. ( 4.21 ) becomes independent of x: 



kl 



(4.30) 



where the constant was evaluated as C^r —i{2T:)'^{NZ^)^^/^ . On the other hand, Eq. ( 4.22 ) shows a narrow peak 
at X = 1/2. Therefore, the pseudoscalar state is a highly collective state, while the scalar state shows an approximate 

constituent picture. 

Now let us compare our result Eq. (4.21) with those of the literatures [ p^p] ]. First of all, equiv alence with the 
result of Ref. [Tl| ] is easily verified. As we commented before, the unfamiliar spinor structure in Eq. (4.16) is due to 
our specific choice of the mode expansion and the representation for the 7 matrices. If one uses the following mode 
expansion for the good component of the fermion: 



2tt 



dk^ 



V27rfc+ 



6(fe, A)u+(fe, X)e'^^ + d^k, X)v+{k, A)e 



-ikx 



one can obtain the same spinor structure as that of Ref. ||Tl|]. Of course the two LC wavefunctions should 
coincide with each other for observable quantities. Indeed, both give the same (quark) distribution function 

q(x)=/d2fc^/(2^)3E„_^|$„^(fc)p . 

On the other hand, the result of Ref. jlOj seems different from ours Eq. (4.21). The possible origin of the discrepancy 
might be attributed to the following two points. First of all, the author of Ref. ||l^ considered the Melosh transforma- 
tion ||2^ which relates the LF spinor and the usual spinor in the equal-time quantization. Such nonstatic spin effects 
might be important when we discuss phenomenological aspects of light mesons (for example, see Ref. [p3[). However, 
even if we take it into account, it is hard to see the coincidence. Secondary but most importantly, he derived the pion 
LC wavefunction by projecting the Bethe-Salpeter amplitude on the equal LC time plane. Though this procedure 
should give the same result as that of the LF bound-state equation as far as we are considering only the ladder (1/A^ 
leading) contribution, equivalence of the two is a highly nontrivial problem in our complicated analysis. 



C. The Gell-Mann-Oakes-Renner and PCAC relations 



Now that we have the LC wavefunction for the pseudoscalar meson, it is straightforward to obtain the decay 
constant f.^: 

^P''U^{0\J^{0)\7^■,P). (4.31) 

For actual calculation, it is safer and easier to treat the plus component. If we use the extended parity-invariant 
cutoff, the result is 
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(4.32) 



Together with the pseudoscalar mass ( 4.26 ) in the chiral hmit, wc find the Gell-Mann, Oakes, and Renner relation, 

NM\ 



9o 



-AniQ (0| |0) 



(4.33) 



The PCAC relation is also checked by using the state | 7r;P). If we normalize the pseudoscalar field TTn{x) cx 
^'(x)«75\['(x) as 



(0| ^n(O) I vr;P) = 1 



— 1/2 /i 

we find that Zt^ given in Eq. (4.27) is the normalization factor 

where we have used the gap equation. Therefore, we arrive at the PCAC relation 

df,j^ = 2mo^'(x)i75*(a;) 
= mlUTTnix) . 



(4.34) 



(4.35) 



(4.36) 



Note that the decay constant ( 4.32 ) and the normalization factor ( 4.35 ) are equivalent to the previous results 
(Eqs. (5.25) and (5.28) in Rcf. M) in the infinitely heavy mass limit of bosons ^ ^ oo. 



D. Symmetric phase 

Her e we consider the symmetric phase in the chiral limit mg = 0. When < g^^ — 27r^/A^, the gap equation 
( |3.26 ) has only a trivial solution M — 0. A quantity which should be zero in the broken phase is now estimated as 



Subsequently the factor a defined by Eq. ( D13 ) is different from that of the broken phase [Eq. ( 4.20 )], 

-1 _ -1 1 9o , 2gg dx [ 2 qo^ 

Then, both of the LF bound-state equation for the scalar and pseudoscalar states are given by 

<.^..(.,fcl) = ^0s..(.,fcl)-f^^f / d^i.^^KM) . (4.39) 
The solution to the bound-state equation is 



2 / „2 



where Cgym is a normalization constant and rrisym = Ti^r — is given as a solution of the equation 



1 1 r' r d'k 



72/,. Tr)2 



- —= dx f4 411 

gl gl io ""J {2^fkl-m%^A^-x) ' > 

Again this is equal to the previous result of the chiral Yukawa model with /x^ — > 00 (Eq. (5.26) in Ref. [Q) and therefore 
if we use the same cutoff as before, we obtain the same result for msym- Moreover, though the above calculation was 
intended only to g^ < g^^. case, if we increase the coupling constant over its critical value g^^, we find a negative 
solution TOgyjjj < 0. This implies that the symmetric solution causes instability when g^ > g^^. and thus we must 
choose the broken solution. 
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V. SUMMARY AND CONCLUSION 



We have investigated a description of D^SB on the LF in the NJL model. The importance of solving the fermionic 
constraint for the bad spinor component was greatly stressed in analogy with the zero-mode constraint of scalar 
models. The exact classical solution enabled us to check the properties of the LF chiral transformation. Though 
the chiral transformation is differently introduced on the LF, we finally found the equivalence to the usual chiral 
transformation. 

For a description of D^SB of LF NJL model, it was very important to solve the fermionic constraint nonperturba- 
tively in quantum treatment. To do so, we introduced a bilocal formulation and solved the bilocal fermionic constraint 
with a fixed operator ordering by the 1/iV expansion. Systematic expansion of the fermion bilocal operator is 
realized by the boson expansion method as the Holstein-Primakoff expansion. The leading bilocal fermionic constraint 
was identified with the gap equation for the chiral condensate after we took care of the infrared divergence. If we 
choose a nontrivial solution of the gap equation, we have a Hamiltonian in the broken phase but with a trivial vacuum. 

The physical role of the fermionic constraint in the LF NJL model is very similar to that of the zero-mode constraint 
for scalar models. We have seen a close parallel between these two constraints. Especially it should be noted that the 
gap equation came from the longitudinal zero mode of the bilocal fermionic constraint. 

It is very natural that we can reach the broken phase by solving the quantum fermionic constraint by 1 /N expansion 
because the fermionic constraint is originally a part of the Euler-Lagrange equation and thus must include relevant 
information of dynamics. What we did is very similar to the usual mean- field approximation for the Euler-Lagrange 
equations. Indeed the leading order in the expansion corresponds to the mean-field approximation. However, our 
way of solving the fermionic constraint with the boson expansion method can easily go beyond the mean-field level. 
Such higher order calculation enabled us to derive a correct broken Hamiltonian and to show the divergent behavior 
of the (spatial integration of) pseudoscalar field. 

Independence of mass from the mode expansion has both desirable and undesirable aspects. The inclusion of correct 
mass dependence into the IR divergent integral was required when we identify the lowest fermionic constraint with 
the gap equation. This is the point we must always take care of. On the other hand, the Fock vacuum is defined 
independent of the value of mass. Due to this fact, it is enough to have only one vacuum, namely, the Fock vacuum 
even in the chirally broken phase. This is the favorable aspect. However, the cost of such a simple vacuum was payed 
by, for example, unusual chiral transformation of fields such as [Q\^ ^ ^'175^'] 7^ —21^"^ and non-vanishing of the LF 
chiral charge \Q\^ ,H] 7^ in the broken phase. We found that the both effects are proportional to the dynamical 
fermion mass M. We also insisted the necessity of a bare mass term which accurately produced the constituent mass 
term. Although the special role of the fermionic constraint might be restricted to the LF NJL model, the unusual 
chiral transformation and the nonconservation of the chiral charge are general features of the chiral symmetry breaking 
on the LF. This is because they are natural consequences of the coexistence of the chiral symmetry breaking and the 
Fock vacuum. 

The leading order eigenvalue equation for a single bosonic state is equivalent to the leading order fermion-antifermion 
bound-state equation. The bound-state equations were solved analytically for scalar and pseudoscalar mesons and 
we obtained their light-cone wavefunctions and masses. The meson masses, the decay constant, and so on were fairly 
consistent with those of our previous analysis on the chiral Yukawa model. The leading order calculation was limited 
only to two body sector (fermion and antifermion) . If we consider the higher order Hamiltonian such as /i*-^-* or h^'^\ 
we will be able to discuss four or six body sectors. In other words, since we have bosonic meson states, we can expand 
the Fock space in terms of the mesonic degrees of freedom. Then, for example, we will be able to discuss the mixing 
of scalar state and two pseudoscalar fields (tt-tt mixing with a). 
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APPENDIX A: CONVENTIONS 

We follow the Kogut-Soper convention |Q . First of all, the light-front coordinates are defined as 
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x^^^{x°±x^), x\=x' (i = l,2), (Al) 
V 2 

where we treat x'^ as "time" . The spatial coordinates x~ and x± are called the longitudinal and transverse directions 
respectively. Derivatives in terms of are defined by d± = d/dx"^. It is useful to introduce projection operators K± 
defined by 

A± = = ;7|7"7^. (A2) 

Indeed A± satisfy the projection properties Aj. — A±, A_|_ + A_ = 1, etc. Splitting the fcrmion field by the projectors 
as 

^^i" = %l;l + iP'L, ipl = Ai^", (A3) 

we find that for any fermion on the LF, ?/^_ component is a dependent degree of freedom. -0+ and "0- are called the 
"good component" and the "bad component" , respectively. 

As is commented in the text, for practical calculation, we use the two-component representation for the gamma 
matrices. The two-component representation is characterized by a specific form of the projectors (1.2). Then the 



projected fermions ^p± have only two components. There are many possibilities which realize Eq. (1.2). For example, 
a specific representation 



is used in Ref. [2l] . In this paper, however, we choose a representation ( |l.4| ) from which it is easy to extract information 



of the 1-1-1 dimensional results. Two-component spinors V' and x are defined by Eq. (1.3). Results of the chiral Gross- 
Neveu model can be easily obtained if we make a replacement for the Pauli matrices (73 —s- 1 and cr* — s- 0, and regard 
-0 and X as one component spinors. 

Using this representation, the Lagrangian density of the NJL model is written as 

£ = zij^d+i: + ix^d-x - ^ (V-V^a.x - x^^^s,^) - ^ {^h + xH) 

+ Y { (^^x + X^V-)' - (V-V^x - xVV)2} . (A5) 



APPENDIX B: CLASSICAL SOLUTIONS TO THE FERMIONIC CONSTRAINTS 

To solve the fermionic constraints classically means that we treat all the fermion fields (both good and bad 
components) as Grassmann numbers and neglect all the c-numbers which will emerge in quantum theory under the 
exchange of variables. 

Before discussing a complicated equation of the NJL model, it would be better to go first with the chiral Gross- 
Neveu model. We solve the fermionic constraint of the chiral Gross-Neveu model with the antiperiodic boundary 
condition: 

[id^+g\{x-)]x=^i^ , (Bl) 

Xa(x^ = -00) = -Xa{x^ = 00) , (B2) 

where we used a matrix notation with a matrix aab{x~) = 4'a{x~)ipl{x~). The solution to this equation is given by 

/oo 
dy^GcNix^.y^)-^ i'iy^) , (B3) 

where Ggn{x^ ,y ) is the Green function satisfying 

{tdl+g\{x-)} GGN{x-,y-) = S{x- - y-) , (B4) 
Ggn{x^ = -00, y") = -Ggn{x^ = 00, 2/") . (B5) 
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Due to Eq. (B5), the solution of course satisfies the antiperiodic boundary condition. Equation (B4) is solved as 



GGN(a; ,y ) = g''qI{x ) 



C 



G 



(0)- 

GN 



\y-) 



Gg)(x-) = Pe-- 



C = - 



1 (ex.) 



G, 



(0) / 



(B6) 

(B7) 
(B8) 



-00 ) 



where GqI^{x^) is a solution of a homogeneous equation {id^ + g^a{x^)] GqI^{x^) ~ and the integral constant G 
has been determined so that Gq,^[x^ ,y~) satisfies the antiperiodic boundary condition. When = 1, the solution 
(B3) is equivalent to Domokos' solution to the Thirring model on the light front p^ . 

In two dimensions, the LF chiral transformation is not distinguishable with the C/(l) transformation. Indeed, the 
"LF chiral transformation" on the good component is t/i ^ 6*^-0 f^nd equivalent to the ?7(1) transformation. And also, 
the solution (B3) implies that the bad component rotates just the same way as the good component x 
Next let us turn to the NJL model. If we explicitly write all the indices, the fermionic constraint (2^1 



■ e 
is 



id- 



X2a 



Ipla^lt Xlb - 01a'02b Xlb 
V'2qV'26 X2b - 1p2a1plb xlb 



(B9) 



where dz = di — 182 and dg = di + 182- Since the equation for xi (or X2) includes xi and (o^' X2 and x|), it is 
useful to introduce a constraint equation for —X2 instead of X2- Then we have a more tractable equation 



Xla 

-Xla 



1 



m^lpla - dz^2a 



v^LV'2& 



Xlb 

A26 



As in the 1+1 dimensional case, the solution is immediately given if we find the Green function G(x , y 
satisfies 



(BIO) 
, x± ) which 



with a matrix Aija bjx) defined by Eq. (| 
Eqs. (O) and (O) in the text. 



{idl + g^Aix-)}G{x- ,y- ,x^_) ^ 5{x- ^ y-) , (Bll) 
). The result is very similar to the two-dimensional result and is given by 



APPENDIX C: PROBLEM OF THE OPERATOR ORDERING 



Here we consider the problem of operator ordering within the chiral Gross-Neveu model with = 1. Following the 
standard procedure, the Dirac brackets are calculated as 



{i^{x),ipHy)}^ = -iS{x -y ) 
{x(x),ip^{y)}^ -iGGN{x,y) 



{x{x),Tp{y)}^ 



mo _ 
'iGQ^{x,y)g'^i}X{y) 



9^4{y)x{y) 



(CI) 
(C2) 
(C3) 



where Ggn(2:,?/) is the Green function (B6) for A^ = 1 case. To quantize the system we simply replace the Dirac 
bracket {A,B}y) by the anticommutation relation —i{A,B}. This procedure has ambiguity of the operator ordering. 

The operator ordering we took for the fermionic constraint (2.1) in the NJL model corresponds to the following one 
in the chiral Gross-Neveu model. 



id-x + 9 WX 



t _ mo 



V2 



We can easily find its quantum solution due to [ipip'^{x)^ ipip'^ (y)] — 0. The solution is 

/oo 
dy~ GGNix,y)—^ilj{y) , 



(C4) 



(C5) 
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where G is again the Green function (B6) with N = 1. 

Now let us consider the consistency for the anticommutator {XiV"^}- It can be cal culated two different ways: (i) 
from the solution Xsoi of the fermionic constraint, and (ii) from the Dirac bracket (^2|). We fix the operator ordering 
of the fermionic constraint by Eq. (C4) and check whether the Dirac bracket can produce the same anticommutator 
or not. 

Instead of the anticommutator itself, we present here the calculation of a quantity iD'-^{x{x) , where iD'-^ = 

id- + g'^tpifj^ . Using the solution (C5), we have 



iD^l {Xso^ix), i^Hy)} = 6{x- - y-) ~ gV^x) ■ 



(C6) 



On the other hand, if we take the simplest ordering in the r.h.s. of the Dirac bracket (C2), we obtain 



iD't{x{x),i^\y)] =iDlGGN{x,y) 



Too 
V2 



Too 



S{x -2/ ) ( ^ -5V'^x(y) 



(C7) 



This is identical with the result (C6). Therefore we find our ordering Eq. ( |C4[ ) is consistent with the anticommutation 
relation 



{xix),tpHy)} = GGNix,y) ( ^ - ffV^X 



(C8) 



Of course if we take other operator ordering, the two results do not coincide. We expect that even in the NJL model, 
we c an s elect the r.h.s. of Dirac brackets so that they coincide with the direct result with the ordering defined by 
Eq. (O). 



APPENDIX D: BILOCAL FERMIONIC CONSTRAINTS AND THEIR SOLUTIONS BY THE BOSON 

EXPANSION METHOD 



It is tractable to solve the equations for T{x,y) aiidU{x,y) rather than 7^^(a;,y) and ^-/^^(a;, y). In momentum 
representation, the fermionic constraints for T and U are 

1 



(Dl) 



f d^p'd^q' 



q+i U{p,q) = - 



2 J (2^)3 {^o.0ip,q-p'-q'){Sapr{p',q') + al^iUip',q')) 

- (So^pTip', q') - cj%i u(p', q')) MMq p' q', p)} , 

{<^^i-iqW +ma)}^^Ma0{p,q) + {{-iq'j^cr' + mQ)a^} ^^Mai3{q,p) 
'2 J ~(2^ 



{Mapip, q p' q') K^rip', q') + 6^pi U{p' , q')) 



{alpTip', q') - 5^pi U{p', q')) M^p{q - p' - q' , p)} 



(D2) 



In place of the quantization condition ( |3.l| ), the system with bilocal operators can be characterized by the following 
algebra 

■ Map{Pi,P2) ■■ , : -^75(91.92) 

= : Mas{Pi,q2) ■ Sp^S'^^Hp2 + qi) - ■ M^f3{qi,P2) ■ S^sS^^Xpi +92) 
+N5^s5pi5^'^\Pi+(l2)5^^\P2 + qi) 

X {e{pt)0{pt)0{-ql)e{-q+) - 0{-pt)9{-pt)e{qt)9iqt)) , (D3) 

where the normal order of M was defined with respect to the Fock vacuum (|3.3|) 
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■.M+7{p,q): |0) = :A4-+(p,q): |0) = : q) : |0) - 



(D4) 



The upper indices stand for signs of the fongitudinal momenta. 

The comphcated structure of the algebra for the bilocal operators which originates from the fermion statistics, is 
greatly reduced if one introduces the boson expansion method. We can represent the operators : : in terms of 
bilocal boson operators B{p, q) of order 0{N^) so that they fulfill the original algebra (D3). Since the algebra has a 
bosonic feature in the large iV limit, 



it would be better to choose a representation which satisfies this|. The Holstein-Primakoff type expansion satisfies 
the requirement. 

Physically this procedure corresponds to extracting purely bosonic degrees of freedom from a fermion-antifermion 
system, i.e., a mesonic system. The power of the boson expansion method in the light-front field theories was first 
demonstrated by one of the authors He applied the Holstein-Primakoff type expansion to 1-1-1 dimensional QCD 
and derived an effective Hamiltonian for mesons as local bosons whose masses are given by the 't Hooft equation. 
Using the effective Hamiltonian, we can in principle study, say, scattering of mesons as qq bound states. 

Since the essential structure of the algebra (D3) is determined only by the longitudinal momentum, it is straightfor- 
ward to apply the Holstein-Primakoff expansion discussed in Ref. ]l5[ | to four dimensional fermionic theories. Indeed 
the operators : M : are represented as follows: 



■Ml+{p^,P2) 
■M-^{pi,P2) 



/ [dq]J2^U(^Pi^l'>^l^'riP2,Q) = -4^a(P2,-Pi) 

- / M9]Xl-^7/3(9>-P2)57«(g>Pl) > 

7 

/ [dq] ^(ViV - A)fi^{p2, q)B^a.{q.Pi) , 
•' 7 

I [dq]Y,B\p{q,-P2){VN^)^„{q,-p,) . 



(D5) 
(D6) 
(D7) 
(D8) 



These give the 1/7V expansion of Aiapip, q)- The first few terms are shown in the text [Eqs. (|3T0|)-(|3l^)]. 

If we expand the bilocal operators T{p,q), U{p,q), and Maf3{p,q), the equation for order n can be written in a 
compact form: 



u(^\p,q) lG(")(p,g) 



2f(pl) 

9o „ I 



d^k ( t^^){k,p + q-k) 
(27r)3 U("Hfc-P + 9-fc) 



(D9) 



where quantities j^^"^ (p, q ) and G'^'^^{p,q) are generally complicated functions of bilocal operators except for the 
lowest order [see Eq. ( |3.18[ )]. For example, F^^'' and G*^^^ are 



F^'\p,q) = T^^HqW + MU UlliP.d) &Q,P)) , 



GW(p,q) 



2q+ 
—i 
2^ 



(1) 



(DIO) 
(Dll) 



where fi''^^ {p, q) is given by the boson expansion method Eq. ( [3.11 ) . Since all of the orders of the operators are less 
than n, we can solve this equation order by order. The solution of this integral equation is 



t(")(p,q) 
q) 



F^"Hp,q) 
G^^\p,q) 



<Pk' f F(-^'){k\p + q-k') 



T?^"^^^ ^'^^L^WrK G(") {k\p + q-k') 



(D12) 



^Actually there are many possibilities to express Eq. (^^) in terms of bosonic operators, corresponding to various "local 
expansions" of the Grassmannian manifold of the bilocal operators. 
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where 



(2^)3 P+ - k+ 

The quantities t^'^^{p, q) and u'-^\p, q) are necessary for obtaining a correct Haniiltonian of the system. 



-iP^)-{^+9ll 7S3#^) • (D13) 
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